The  besy1 function returns the y1st Bessel function of its argument.
  besy1 expects its argument to be in radians.